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Abstract 

We  present  a  phenomenological  model  that  indicates  with  high  probability  that  the 
smaller  of  the  two  energy  gaps  in  MgB2,  the  so-called  n  gap,  contains  nodal  lines  with  a  six-fold 
symmetry  (i-wave).  The  model  also  indicates  that  the  larger  gap,  the  so-called  a  gap,  is  a 
conventional  s  wave.  The  model  is  an  extension  of  the  BCS  theory  that  accounts  for  the  elastic 
anisotropy  in  MgB2  and  the  Coulomb  repulsion.  It  is  based  on  a  phononic  pairing  mechanism 
and  assumes  no  coupling  between  the  two  energy  gaps  in  MgB2  at  zero  temperature.  All  of  the 
parameters  of  the  model,  such  as  the  sound  velocities,  are  independently  detennined  material 
constants.  The  results  agree  with  a  previous  ad-hoc  hypothesis  that  the  n  energy  gap  has  six 
nodal  lines.  That  hypothesis  was  motivated  by  low-temperature  measurements  of  the  surface 
impedance  and  intermodulation  distortion  in  high-quality  thin  films.  We  briefly  review 
experimental  evidence  in  support  of  our  hypothesis  and  review  evidence  in  the  literature  that  has 
led  to  the  conclusion  that  both  gaps  are  s  wave.  We  find  that  the  evidence  from  the  literature  for  s 
wave  is  inconclusive.  Our  finding  is  that  the  n  gap  has  six  nodal  lines. 
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1  Introduction 

Superconducting  magnesium  diboride  has  been  demonstrated  to  possess  properties  that 
are  promising  for  applications  at  microwave  frequencies.  The  relatively  high  transition 
temperature  of  40  K,  low  surface  resistance,  relatively  long  coherence  length,  and  the  high 
critical  magnetic  field  of  >1  T  have  led  to  interest  for  microwave  electronic  applications,  such  as 
filters  and  for  superconducting  radio  frequency  (SRF)  applications  in  accelerator  cavities[l]. 
MgB2  has  also  attracted  fundamental  interest  because  of  the  two  superconducting  energy  gaps, 
while  at  the  same  time  it  appears  to  be  a  superconductor  where  the  pairing  mechanism  is 
phononic  as  in  BCS  theory  [1,2]. 

In  investigations  of  the  microwave-frequency  properties  of  high-quality  thin  films  of 
MgB2,  we  have  previously  found  strong  evidence  that  the  smaller  of  the  energy  gaps  has  a 
symmetry  that  includes  nodal  lines,  and  the  symmetry  is  six  fold  (£  =  6,  i-wave  in  spectroscopic 
notation)  due  to  the  hexagonal  crystal  structure  of  MgB2, 

A<r(^7’)  =  A<r(r)sin(6p)»  (1) 

where  An(T)  is  the  n  energy-gap  magnitude  and  (p  denotes  the  azimuthal  angle  in  the  cylindrical- 
coordinate  system  of  the  relative  k- space.  This  proposition  has  been  reported  in  series  of 
papers[3,4,5].  These  findings  will  be  discussed  in  more  detail  in  Section  II.  Recently, 
measurements  of  the  microwave  properties  of  MgB2  by  another  group  [6],  using  a  different 
method  and  with  different  high-quality  films,  have  been  reported.  As  the  analysis  in  Section  II 
demonstrates,  this  experimental  evidence  agrees  with  the  i  =  6  hypothesis,  while  inconsistent 
with  s-wave  symmetry. 

To  give  the  t  =  6  hypothesis  a  theoretical  underpinning,  we  have  developed  a 
phenomenological  model  for  calculating  the  energy  gaps  in  MgB2  at  7=  0  K.  The  presentation  of 
this  model  is  the  main  purpose  of  this  work.  The  results  strongly  support  the  hypothesis  that  the 
smaller  gap,  the  so-called  n  gap,  exhibits  i-6  symmetry  while  the  larger  gap,  the  a  gap,  has 
conventional  s-wave  symmetry.  The  model  is  parameter-free,  meaning  that  all  parameters  used 
are  available  from  other  calculations  or  measurements  in  the  literature.  The  model  is  inspired  by 
previous  work  regarding  the  d-wave  symmetry  of  the  energy-gap  of  YBCO,  a  cuprate  [7],  This 
work  builds  on  a  previous  generalization  of  the  BCS  theory,  where  the  coulomb  repulsion  is 
included  [8],  and  adds  the  important  feature  of  large  anisotropy  in  the  elastic  properties.  Our 
work  differs  from  [7],  however,  in  several  respects.  It  is  applied  to  MgB2  where  the  pairing 
mechanism  is  believed  to  be  phononic  [9,10,11],  the  attractive  interaction  parametrization  is 
distinct,  and  it  aims  at  the  f=0K  energy  gap  rather  than  a  focus  on  the  T~TC  domain 

We  are  aware  that  it  has  been  widely  concluded  in  the  literature  that  both  energy  gaps  of 
MgB2  have  a  conventional  s-wave  symmetry.  In  the  discussion,  Section  IV,  we  examine 
literature  for  probes  that  imply  the  n  energy-gap  symmetry.  Our  finding  is  that  the  pertaining 
literature  that  suggests  the  s-wave  symmetry  is  inconclusive  and  hence  the  proposition  of  a  i  =  6 
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symmetry  for  n  energy  gap  is  not  inconsistent  with  the  literature.  In  that  Section  it  is  also 
speculated  on  other  instances  where  the  i-wave  energy  gap  symmetry  may  be  realized  and 
propose  further  experimental  work  to  settle  the  symmetry  proposition  of  this  work. 

2  Review  of  the  evidence  for  t-6 

2.1  Intermodulation  Distortion 

As  reported  previously  [3],  microwave  measurements  of  the  temperature  dependence  of 
the  intermodulation  distortion  (IMD)  in  MgB2  at  low  temperatures  show  a  1/J2  increase  as  the 
temperature  is  lowered.  As  a  reminder,  Figure  1  shows  this  previously  reported  result,  where  the 
similar  data  from  YBCO  is  included  [12]. 

The  1/72  behavior  in  YBCO  is  attributed  to  the  nonlinear  Meissner  effect  [13,14].  It  is  a 
signature  of  a  nodal  energy  gap,  a  d-wave  in  the  case  of  YBCO.  In  analogy  to  YBCO  case,  we 
have  interpreted  the  observed  low-temperature  increase  in  IMD  in  MgB2  as  a  strong  indication 
that  the  energy  gap  contains  nodal  lines.  As  discussed  in  detail  in  Section  III,  the  hexagonal 
crystal  symmetry  of  MgB2  and  the  observed  \/T2  behavior  of  IMD(7)  implies  the  i-  6 
symmetry.  For  an  s-wave  energy-gap  symmetry  the  IMD  vanishes  exponentially  as  T— >0  K.  No 
other  plausible  explanation  has  been  offered  for  the  measured  IMD(T). 

2.2  Penetration  Depth  A A,(T) 

Another  piece  of  evidence  for  the  t—  6  symmetry  is  the  low-temperature  variation  of  the 
penetration  depth  reported  and  analyzed  in  Ref.  [4].  Below  the  temperature  of  about  5  K  the 
penetration  depth  increases  as  the  temperature  is  lowered.  Figure  2  shows  the  result  reported  in 
Ref  [4]  of  fits  that  assume  i  —  6  symmetry  for  the  dominating  n  energy  gap.  We  have  interpreted 
the  increase  as  a  manifestation  of  Andreev  surface-attached  states  that  are  associated  with  the 
existence  of  nodal  lines  in  the  energy  gap.  An  identical  observation  was  reported  for  YBCO, 
which  has  nodal  d-wave  energy-gap  symmetry  with  two  nodal  lines  [15,  16]. 

2.3  Surface  Resistance  vs  Temperature  RS(T) 

The  temperature  dependence  of  the  surface  resistance  RS(T)  also  sheds  light  on  the 
symmetry  of  the  energy  gap.  As  discussed  in  the  following,  for  s-wave  symmetry  RS(T)  vanishes 
exponentially  as  the  temperature  is  lowered,  while  for  a  nodal  energy  gap  it  follows  a  power  law. 
Our  data  shows  indeed  a  power-law  temperature  variation  at  low  temperatures,  consistent  with 
the  t  =  6  premise  [5].  Recently,  measurements  [6]  of  the  RS(T)  for  high-quality  MgB2  films  from 
a  different  source  from  that  of  [3,4]  and  by  a  different  measurement  method  have  been  published 
by  a  group  at  Jefferson  National  Accelerator  Facility  (JLab).  As  demonstrated  in  the  following, 
our  new  analysis  of  these  measurements  confirms  that  they  too  support  the  i  —  6  hypothesis. 

The  presumed  symmetry  of  the  MgB2  a  energy  gap  as  an  s-wave  [1]  implies  that  at 
sufficiently  low  temperatures  the  n  energy  gap  dominates  the  temperature  dependence.  Hence,  in 
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the  low-temperature  domain  and  in  the  clean  limit,  an  s-wave  n  energy  gap  yields  a  surface- 
resistance  R s  that  obeys  the  BCS-theory  textbook  expression  that  is  characterized  by  an 
exponential  temperature  dependence,  specifically, [17, 18] 


RS{T)  =  RRES+Cf2pe-A(T)P .  (2) 

In  Eq.  (2)  /  is  the  frequency,  P  =  l/(kB7)  where  kB  is  the  Boltzmann  constant  and  Rres,  C  are 
temperature-independent  fitting  parameters  of  the  residual  surface  resistance  and  a 
proportionality  factor,  respectively.  The  quantity  Rres  is  experimentally  extracted  from  the  very- 
low  temperature  surface-resistance  data,  where  no  measurable  temperature  dependence  is 
observed.  The  origin  of  this  tenn  is  presumed  not  related  to  bulk  superconductivity  [19].  Details 
of  the  JLab  experiments  have  been  given  elsewhere  [6].  In  brief,  the  measured  single-crystal 
MgB2  films  were  grown  by  the  HPCVD  method  [1],  and  have  been  measured  in  a  TEon 
cylindrical  Nb  cavity  with  a  sapphire  rod  inside.  The  resonant  frequency  is  7.4  GHz.  The 
analysis  of  the  loss  and  shift  in  resonance  frequency  in  this  setup  yields  the  effective  surface- 
resistance  and  penetration-depth  dependencies  on  temperature. 

The  Rs(T)  expression  in  Eq.  (2)  suggests  a  semi-log  plot  vs.l/T.  In  such  a  plot,  data  that 
conforms  to  Eq.(2)  would  lie  along  a  straight  line  with  slope  proportional  to 
AK(T  =  0K)»  2.3  meV,  the  accepted  value  for  the  n  energy  gap  [1,2,20].  The  plot  in  Figure  3) 
shows  for  clarity  only  the  results  for  the  350-nm-thick  film  from  [6].  Similar  results  are  obtained 
for  the  200-nm-thick  film.  This  figure  convincingly  shows  that  the  data  does  not  follow  the 
expected  straight  line  for  an  assumed  s-wave  energy  gap  value  AK(f=  0  K)  »  2.3  meV.  As  a 
consistency  check  to  the  analysis  of  the  Rs  data  in  Figure  3),  we  plot  in  Figure  4)  the  Nb  data, 
also  reported  by  the  JLab  group  [6],  As  expected,  the  BCS  expression  fits  the  low-temperature 
data  using  the  literature  value  of  the  Nb  energy  gap,  ANb(r=  0  K)  =  1.5  meV  [21]  for  about  two 
orders  of  magnitude. 

The  inconsistency  of  the  low-temperature  surface-resistance  data  with  the  BCS 
expression,  Eq.(l),  demonstrated  in  Figure  3)  is  resolved  by  analyzing  the  data  under  the 
previously  proposed  assumption  that  the  n  energy  gap  has  a  symmetry  with  six  nodal  lines, 
Eq.(l).  This  energy-gap  symmetry  yields  for  the  surface  resistance  [3] 


RS(T)  =  R 


RES 


P 


15^  1 

8 


V  °  J 


p 


-  =  _A,{T  =  0K) 


kBT 


(3) 


where  a  >  0  and  Rres  are  defined  above.  Since  at  low  temperatures  the  variation  of  AK(T)  is 
negligible  we  approximate  in  Eq.  (3)  An(T)  «  A JT  =  0  K)  [20].  The  expression  in  Eq.  (3)  is  based 
on  the  two-fluid  model,  where  at  low  temperatures  the  temperature  variation  reflects  primarily 
the  corresponding  temperature  variation  of  the  quasiparticle  density  [21].  Note  the  qualitatively 
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distinct  temperature  variations  of  expressions  Eq.  (2)  and  Eq.  (3),  an  exponential  vs.  a  power- 
law. 


Figure  5)  shows  the  results  of  fitting  the  Rs  (T)  -  Rres  of  the  350-nm-thick  film  from  JLab 
with  the  expression  in  Eq.(3).  For  comparison,  measurements  made  with  stripline  resonators  at 
MIT  Lincoln  Laboratory  are  also  plotted  [22],  In  order  to  plot  the  stripline  data  at/=  2.2  GHz 
and  the  Jlab  data  at/=  7.4  GHz  on  the  same  graph,  the  MIT  Lincoln  Laboratory  measurement 
was  scaled  by  the  usual  / 2  dependence.  The  Lincoln  Laboratory  measurements  are  of  a  film 
deposited  by  reactive  evaporation  [23]  and  unpublished  measurements  of  a  film  deposited  by 
HPCVD  process  which  is  also  the  case  for  the  films  measured  at  JLab.  As  Figure  5) 
demonstrates,  the  fits  are  very  good  in  all  cases.  The  fits,  however,  show  variability  in  slope  and 
fitting  parameters.  Such  variability  has  its  probable  origin  in  uncontrolled  factors  in  the  film 
deposition  method,  and  is  normal  in  measurements  of  superconducting  thin  films.  The  point  of 
Figure  5)  is  that  regardless  of  such  extraneous  uncontrolled  factors,  all  films  share  the  same 
particular  temperature  variation  of  Eq.  (3).  This  outcome  resolves  the  inconsistency  encountered 
in  Figure  3).  We  emphasize  that  the  Jlab  measurements  and  above  analysis  are  not  a  stand-alone 
proof  of  the  unconventional  energy-gap  symmetry,  Eq.  (1),  but  rather  provide  additional  support 
to  the  original  proposition  that  was  based  on  measurement  of  the  IMD  [3].  Accounting  for  the  a 
energy-gap  contribution  at  higher  temperatures  lies  outside  the  scope  of  this  work. 


3  Phenomenological  model  for  the  energy-gap  symmetry 

In  this  Section  we  show  that  the  £  =  6  energy-gap  symmetry,  Eq.  (1),  is  also  substantiated 
by  calculations  based  on  a  phenomenological  model.  That  model  was  introduced  for  and 
successfully  applied  to  YBCO  [7],  where  d-wave  symmetry  of  the  energy  gap  is  supported  by  a 
large  body  of  data  [24].  The  model  is  based  on  a  phononic-pairing-mechanism  premise,  which  is 
generally  accepted  to  be  the  case  for  MgB2  [9].  It  generalizes  the  classic  BCS  theory  in  two 
respects.  It  accounts  for  the  mutual  Coulombic  repulsion  between  the  two  electrons  in  a  Cooper 
pair  [8]  and  for  anisotropy  in  the  superconductor  elastic  properties  In  YBCO  [7]  and  MgBo  this 
anisotropy  is  substantial,  see  below.  All  the  required  input  parameters  are  extracted  from 
pertaining  calculations  and  data  in  the  literature. 

A  complication  in  MgB2  that  is  absent  in  YBCO  is  the  existence  of  two  energy  gaps,  the 
n  (small)  and  a  (large)  energy  gaps,  respectively.  These  are  weakly  coupled  due  to  their  distinct 
underlying  bonding.  At  temperatures  T  ~  Tc  however,  where  the  energy  gaps  become  small  [25, 
26,  27],  this  coupling  cannot  be  ignored.  By  contrast  in  the  T— >  0  K  limit,  it  can  be  neglected 
given  the  considerable  disparity  in  the  energy  gap  magnitudes 
(T  =  0)  «  0.002  eV,  ACT  (T  =  0)  «  0.007  eV  .  Consequently,  at  T=  0  K  the  two  energy  gaps 
can  be  treated  separately,  a  feature  that  simplifies  the  analysis  considerably. 
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Our  objective  is  solving  the  general  BCS  energy-gap  equation  at  7=  OK, 


V 


A  k' 


k' 


,  ,  ,  .  f  dk'v(k,k')  .  . 

2  E[k')  (2<t)!  >  y  2E[k  j 


(4) 


The  symbols  in  Eq.(4)  denote  the  following.  The  energy  gap  at  7=0  K  is  A  |  k  j ,  the  momentum 

is  k  and  the  summation  (integration)  is  over  the  first  Brillouin  zone,  denoted  by  BZ.  The  sample 
volume  is  V  ,  the  chemical  potential  is  p,  and  the  band-energy  and  quasiparticle  energy  are 

is  E^k^,  respectively.  The  mutual  interaction  of  the  two  electrons  in  a  Cooper  pair  is 

v(jc,k'^j  [28].Eq.(4)  has  always  the  trivial  solution,  A^k)  =  0,  which  corresponds  to  a  state 

without  quasiparticles  [28].  We  are  interested,  however,  only  in  nontrivial  solutions  where 
A^kj^O,  if  they  exist.  Technical  details  of  solving  Eq.(4)  are  deferred  to  the  Appendices.  In 
this  Section,  we  focus  only  on  key  ingredients  of  the  model  and  results  of  the  calculations.  In  the 
BCS  theory,  where  A^k  j  and  V  |k,k  ’j  are  assumed  to  be  constants,  the  nontrivial  solution  of 
Eq.(4)  can  be  found  in  textbooks  [21,  29]. 


Consider  first  the  interaction  factor  VyC,k'j  in  Eq.(4).  It  consists  of  an  attraction, 

resulting  from  a  single-phonon  exchange,  and  the  corresponding  Coulomb  repulsion.  When 
accounting  for  realistic  phonon  modes,  this  interaction  factor  is  complex  [30].  Here  we  follow 
the  generalized  BCS  theory  by  adopting  the  following  ansatz  [31,7] 


v  (*»*')  =  ^vph  {4)  o[hcoD  -|£(£)|)  e(naD-\^(k^ 

+  -VC  (q)  e(hcQp-\%(k}  j  d[hcop-\^(k'^  j  =  V^,  |  z(p)  ,  \^(k ') 


(5) 


q=k-k\ 


The  symbols  in  Eq.  (5)  denote  the  following.  The  interactions  V PH  (q)  and  Vc  (q)  denote 
the  attractive  interaction  due  to  a  single  phonon  exchange  between  the  two  Cooper-pair  electrons 
and  corresponding  Coulombic  repulsion,  respectively,  co D  and  (QP  denote  the  Debye  and  plasma 
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frequencies,  respectively,  and  0(g)  is  the  dimensionless  step  function  which  equals  1  or  0  for 

positive  or  negative  j  ,  respectively.  The  step  function  in  Eq.  (5)  expresses  the  BCS  premise 

that  only  electrons  within  a  thin  energy  shell  near  the  Fermi  surface  partake  in  pairing.  That 
energy-shell  width  is  the  Debye  energy.  Similarly,  the  maximum  energy  exchanged  due  to  the 
Coulomb  repulsion  is  the  plasma  energy  [8,31] 

Consider  first  the  attractive  interaction  VPH(q)  in  Eq.  (5).  For  an  isotropic 

superconductor  in  the  weak-coupling  regime,  the  single -phonon  attractive  interaction  is  given  by 
[32] 


VPH  =  -y10(tia>D-ti(o), 


r  = 


(  2  A 
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where  y  is  a  constant  amplitude  (dimensionality  [y]  =  s[&  i 3/2  with  £  and  l  denoting  energy  and 

length,  respectively)  and  flCO  is  the  exchanged-phonon  energy.  The  other  symbols  in  Eq(6) 
denote  the  following:  z  is  the  partaking-ion  valency,  vSOUND  is  the  longitudinal  sound  velocity, 
cjtf  is  the  Thomas-Fermi  momentum.  The  no  parameter  is  the  ion  number  density  (dimensionality 
[n0^j  =  f'3),  i.e.,  zeno  is  the  ionic -background  charge  density  that  neutralizes  the  electron-sea 
charge  density.  The  M  parameter  is  the  partaking-ion  mass  and  B  is  the  bulk  modulus 
(dimensionality  [£?]=££  3  ).  The  second  expression  for  y  in  Eq.(6)  is  obtained  by  employing 
the  spherical  Fermi-sea  expressions  qFF  =  (6/rn0e  2)/sF  ,  sF  =h2kF /(im*), 

kF  =  (()n1n  n  /2 ) '  and  ( vSOUND )  ~  =  b/(M  n  0  )  where  m*  denotes  the  effective  electron  mass.  It 

provides  a  connection  of  y  with  measurable  or  calculated  quantities  in  the  literature.  Eq.  (6)  is 
based  on  the  standard  jelium  model  for  metals,  a  system  comprising  a  positive  ionic-charge 
background,  zeno  that  is  compensated  by  the  Fermi-sea  charge-density  of  electrons  released  by 
the  ions.  Note  that  in  the  relative  distance  between  the  two  Cooper-pair  electrons,  the  interaction 
Eq.  (6),  is  a  delta  function. 

The  point  to  note  in  Eq.(6)  is  the  vSOUND  (or  B )  factor,  which  in  general  depends  on  the 
direction  in  momentum  space.  That  directional  dependence  reflects  the  disparity  in  the  crystal 
elastic  constants  [33].  For  a  hexagonal  crystal  structure  in  particular,  the  longitudinal  sound 

velocities  in  the  w  plane  and  in  the  z  direction  are  given  by  [34] 
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where 


C  j  j ,  C  33  are  the  elastic  constants  in  the  standard  notation  (dimensionality 
=  £  i~2  )  and  p  is  the  ion  mass-density  (dimensionality  [/?]  =  £t"£  5,  where  t  denotes 


time).  In  Eq.(7)  and  hereafter  the  symbols  □  and  _L  denote  quantities  pertaining  to  the  basal 
plane  and  z  directions  in  the  cylindrical  coordinate  system,  respectively.  A  large  disparity 
between  Cx  ,  and  C3  3  implies  a  corresponding  significant  ^-dependence  of  the  interaction  V PH  . 

Deferring  details  to  the  Appendices,  the  parametrization  of  this  dependence  is  a  key  ingredient  of 
our  model.  Specifically  we  use 
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where  the  important  anisotropy  parameter  a  is  defined  [7] 
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In  Eq.(8)  /?(_!_), /?(□)  are  corresponding  bulk  moduli.  Note  that  unlike  for  the  s-wave  case,  in  the 


distance  between  the  two  Cooper-pair  electrons  the  interaction  (8)  is  of  finite  range.  As  Eq.(9) 
implies,  for  an  isotropic  superconductor  the  anisotropy  parameter  vanishes. 


The  Vc  ( q )  factor  in  the  interaction  (5)  is  the  Coulomb  repulsion  between  the  two 

Cooper-pair  electrons.  In  keeping  with  the  jelium  model  for  metals,  we  adopt  the  standard 
Thomas-Fermi  expression, 
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where  the  basal-plane  Fermi-momentum  KF(n),  Kp(o )  refer  to  the  n  and  a  energy  gaps, 
respectively.  While  in  principle  there  are  four  sheets  to  the  Fenni  surface  [9,35],  in  our  model  we 
adopt  just  two  representative  Fermi-surface  sheets,  one  for  each  energy  gap. 

Consider  now  the  energy  gap  A(k)  in  Eq.(4).  It  is  important  to  note  that  its  functional 
dependence  is  severely  constrained  by  the  crystal-structure  symmetry.  As  listed  in  Table  1  [36], 
for  the  hexagonal  crystal  structure  the  only  allowable  multipoles  (in  the  cylindrical-coordinates 
basal  angle)  are  1  =  0,3,6.  Consequently,  the  generalized  ansatz  for  the  gap  function  in  the 
notation  of  Table  1  is  [7] 

ji.  (hlhzY)  for r,* 

fr(kz)  =  <\  forr+  gr{kz)  = 

*z/**z  f°rT  3+ 

(ID 

where  (p  denotes  basal-plane  angle  of  A:  in  a  cylindrical  coordinate  system.  As  the  step  functions 
in  (11)  indicate,  the  gap  A('j  (T)  is  associated  with  the  attractive  phonon-exchange  interaction. 

This  quantity  is  identified  with  the  measured  energy  gap.  The  second  high-energy  gap  A 
in  Eq.  (11)  is  associated  with  the  Coulombic  repulsion  in  Eq.  (5).  In  the  classic  BCS  theory 
/'r  (kz)  -g  v((p)  =  1  and  A^1  (7’]  =  0  .  In  its  generalization  that  includes  a  constant  coulomb 

interaction  [8],  we  have  fr  (kz)  =gT(ycp)  =  1  and  A^(r)^0,  while  in  the  present  calculations 

/r  (Az ) ,  g  v{(p)  have  specified  functional  fonns  dictated  by  the  crystal  structure  symmetry  and 

A(’}  (T)  *  0  .  It  is  expected,  and  numerically  verified  below,  that  A;'*  (T)  A1^1  (T)  in  our 

calculations.  Note  that  by  virtue  of  Eq.(l  1)  the  unknowns  in  the  gap  equation  Eq.  (4)  are  reduced 
to  constants. 

Having  specified  the  interaction,  Eqs.  (8)  and  (10),  and  the  allowable  functional 
dependence  of  the  energy  gap  Eq.  (11),  determination  of  the  pertaining  input  parameters  remains. 
These  are  taken  from  the  literature  and  cited  in  Table  2.  A  few  comments  are  in  order  before 
presenting  the  results.  The  first  concerns  the  elastic  constants  of  MgEE.  To  the  best  of  our 
knowledge  they  have  not  been  measured  [37];  however,  there  are  calculated  values  based  on 
known  band-structure  calculations,  where  the  quoted  values  vary  depending  on  the  methods  of 
calculation  used  [38].  Given  this  situation,  we  adopted  for  our  calculations  average  values, 
specifically 


1  forT  + 

sin(6^)  forTl 

|sin(3^),cos(3^)|  forT 3+ 
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ici,i>ci,2>ci, 3^3,3}  =  {444.3,  60.1,  40.3,  240.8}GPa. 


(12) 


Employing  these  elastic  constants,  Eqs.(7)  and  (9)  yield  for  the  anisotropy  parameter  a  =  0.84, 
comparable  to  that  in  YBCO  [7], 

The  other  comments  pertain  to  the  simplified  Brillouin  zone  and  ionic  parameters  in 
Eq.(8).  Underlying  our  approximation  of  no  coupling  between  the  n  and  a  energy  gaps  is  the 
premise  that  the  phonons  responsible  for  these  two  energy  gaps  are  distinct  [10],  as  are  the 
corresponding  electron  seas.  Indeed,  the  calculated  electron  density  is  accumulated  either  around 
the  B  or  Mg  atoms  [9],  Band  structure  calculations  indicate  that  a  energy  gap  is  associated  with 
vibration  of  the  B  atoms  in  their  graphitic  planes,  in  the  □  direction,  whereas  the  n  energy-gap  is 
associated  with  vibrations  that  involve  the  Mg  atoms,  in  the  _L  direction  [10].  This  consideration 
implies  the  choice  for  the  participating  ions’  valency  parameter  z  and  their  density  no.  A  choice 
for  the  MgB2  unit  cell  is  a  hexagon,  where  the  Mg  atoms  are  located  at  all  its  comers  and  at  the 
centers  of  the  hexagon’s  top  and  bottom  facets.  In  that  hexagon,  the  B  atoms  lie  in  a  hexagonal 
loop  located  in  a  plane  midway  between  the  hexagon’s  top  and  bottom  facets.  A  corresponding 

primitive  unit-cell  is  a  sixth  of  that  hexagon,  of  volume  0.25^3  a2aac  where  a„  and  a,  are  the 

crystal  constants  in  the  basal  plane  and  z-direction,  respectively.  In  this  primitive  unit  cell  there 
is  one  B  atom  and  V2  Mg  atom  with  valencies  1  and  2,  respectively.  Since  the  n  energy  gap 
involves  the  Mg  atoms  it  follows  that  zn  -  2 ,  and  since  the  a  energy  gap  involves  only  B  atoms 

we  set  z  a- 1 .  The  ion  density  parameter,  no  determines  the  electron  number  density,  which 

neutralizes  the  background  ions  charge  density.  Since  the  n  energy  gap  involves  the  Mg,  and 
there  is  !4  such  atom  in  the  primitive  unit-cell  we  set  (n  0)  K  =  V2,  and  a  similar  argument  implies 
that  (n  0)  G  =  1 .  These  choices  ignore  the  different  type  of  bonding  associated  with  the  B  and  Mg 
atoms.  Whereas  the  Mg  atom  is  bonded  (to  the  B  atoms)  by  releasing  its  valence  electrons,  the  B 
atom  is  bonded  (to  B  atoms)  by  sharing  its  valency  electrons  with  its  near  B  neighbor  [39]. 
Notwithstanding  this  consideration,  the  choice  (n  0)  CT  =  1  in  the  calculations  should  be  considered 
as  part  of  our  model. 

The  results  of  the  calculations  are  displayed  in  Figure  6)-8).  Since  our  focus  is  on  the 
£  =  6  symmetry  vs.  the  s-wave  symmetry  of  the  n  energy  gap,  only  these  two  symmetries  are 
considered  here.  Note  that,  in  principle,  combinations  of  different  allowable  symmetries  are  also 
a  possibility.  For  example,  recent  data  suggests  a  small  s-wave  admixture  to  the  dominating  d- 
wave  symmetry  in  YBCO  [40], 

Figure  6)  and  Figure  7)  show  two  calculations  for  the  ratio  of  the  calculated  n  energy  gap 
and  the  experimental  values  as  a  function  the  anisotropy  parameter  a.  The  points  in  the  figures 
are  the  numerically  calculated  values,  while  the  lines  are  a  cubic-spline  fit  to  the  points.  In  Figure 
6)  the  £  =  6  symmetry  is  assumed,  while  in  Figure  7)  the  s-wave  symmetry  is  assumed. 
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Comparison  of  these  two  calculations  for  the  value  of  the  anisotropy  parameter  a  ~  0.83 
suggests  that  the  1-6  symmetry  is  considerably  more  likely  to  be  the  correct  one  than  s  wave. 
This  outcome  is  consistent  with  the  measurements  discussed  in  Section  II.  Once  again  we 
emphasize  that  there  are  no  adjustable  parameters  in  this  calculation.  To  substantiate  our  case, 
the  model  was  also  applied  for  the  a  energy  gap.  We  calculated  this  energy  gap  assuming  both 
the  £  =  6  and  s-wave  symmetry.  These  results  are  shown  in  Figure  8).  For  the  assumed  s-wave 
symmetry  the  calculated  energy  gap  is  within  an  order  of  magnitude  of  the  measured  one.  On  the 
other  hand,  for  an  assumed  t  =  6  symmetry  there  are  no  nontrivial  solutions  to  the  gap  equation 
Eq.  (4)  .  These  outcomes  strongly  indicate  that  the  a  energy  gap  has  s-wave  symmetry, 


consistent  with  current  consensus.  It  was  verified  that  A1,'* (71)  A^(T)  in  all  instances  (not 


d2E 


shown). 


4  Discussion 

Given  that  the  nodal  hypothesis  for  the  n  energy-gap  is  at  odds  with  the  common 
presumption  that  the  symmetry  is  s-wave,  it  is  necessary  to  review  the  literature  that  has  been 
used  to  conclude  and  support  the  s-wave  symmetry  of  the  n  gap  This  has  been  done  in  our 
previous  publications  with  the  conclusion  that,  to  our  knowledge,  the  available  data  does  not 
settle  this  issue  [3,4,5].  We  list  this  review  here,  complemented  by  a  comprehensive  review  of 
microwave  surface-resistance  data. 

The  low-temperature  data  that  was  previously  checked  consisted  of  tunneling,  penetration 
depth,  specific  heat,  muon-spin  rotation,  nuclear-spin  relaxation  (1/7),  Hebel-Schlichter  peak), 
magnetic -field  dependence,  and  the  paramagnetic  Meissner  effect.  In  this  body  of  data  we  find 
that  some  of  it  is  consistent  with  a  nodal  energy  gap  and  some  is  not.  Here  we  compliment  this 
list  by  reviewing  the  available  low-temperature  surface-resistance  data  to  establish  if  it  is 
inconsistent  with  Eq.(3).  There  are  only  a  very  few  previous  reports  in  the  literature  of 
measurements  of  the  microwave-frequency  surface  impedance  of  high-quality  thin  films  of 
MgB2.  As  far  as  we  know  there  are  no  reports  of  single-crystal  measurements.  In  our  opinion  the 
existing  measurements  are  not  definitive  in  distinguishing  the  symmetry  of  the  energy  gap.  The 
following  discusses  the  most  important  of  the  reported  measurements.  A  definitive  result  would 
be  demonstration  of  an  exponential  dependence  of  the  surface  resistance  at  low  temperatures. 

Jin  et  al.  [41],  reported  surface  resistance  measurements  on  relatively  high  quality  thin 
films  and  reported  an  exponential  low-temperature  dependence,  RS(T)M  e  Allykl  as  in  Eq  (2). 

However  the  reported  fitting  parameter  is  A(7=0)  =  3.4  meV,  which  disagrees  with  the  accepted 
values  of  2-2.3  meV  for  the  n  gap,  which  dominates  the  low-temperature  behavior. 
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Zaitsev  et  al.  [42]  reported  a  fit  to  exponential  dependence  over  a  wide  range  of 
temperatures,  but  as  they  point  out  in  the  paper,  the  statistical  errors  in  the  measurements 
preclude  a  definitive  conclusion  of  of  s-wave  behavior. 

Hakim  et  al. [43]  reported  a  fit  to  an  s-wave  exponential  dependence  over  the  entire 
temperature  range.  The  film  was  reported  to  be  a  mixture  of  textured  MgO  and  nanocrystalline 
MgB2,  with  grain  sizes  ~  5  mn  and  T(  =  3  I  K.  While  such  a  film  and  this  result  is  of  interest  in 
the  context  of  a  prediction  that  heavy  doping  of  MgB2  results  in  a  one-gap  high- 7c 
superconductor,  this  data  does  not  pertain  to  the  clean-limit  high-quality  films  under  study. 

Moeckly  et  al.  [44]  reported  measurements  on  a  film  made  by  the  same-reactive 
evaporation  technique  as  we  have  in  our  publications  [3,4,5]  using,  however,  a  parallel  plate 
resonator  as  opposed  to  the  stripline  resonator  that  we  have  employed.  The  results  showed  what 
appeared  to  be  sharp  decrease  at  low  temperature  that  resembled  exponential  behavior  but 
without  further  analysis  to  fit  and  extract  a  value  for  the  energy  gap.  We  point  out  that  the 
parallel-plate-resonator  result  required  subtraction  of  “background  losses”  that  may  have 
increased  the  uncertainty  of  the  measurements.  We  also  note  that  the  paper  also  reports 
measurements  made  in  our  laboratory  on  similar  films  (as  we  have  reported  in  our  papers  and  in 
Section  II)  that  follow  a  power-law  behavior  at  low  temperatures. 

Cifariello  et  al.  [45]  have  measured  the  IMD  in  high-quality  MgB2  films.  Their 
measurements  did  not  show  an  upturn  in  IMD  at  the  temperatures  range  T/Tc  >  0.2.  In  so 
measuring  they  missed  the  region  with  the  strongest  evidence  of  the  upturn  where  T/Tc  <  0.2,  see 
Fig.(l). 

To  conclude  this  survey  of  microwave-frequency  surface-impedance  measurements,  we 
assert  that  the  literature  on  this  topic,  to  the  best  of  our  knowledge,  contains  only  inconclusive 
evidence  for  s-wave  symmetry.lt  is  worthwhile  to  point  out  that  there  have  been  also  low- 
temperature  reports  of  penetration  depth  measurements  that  indicate  unconventional 
superconductivity  [46,  47,  48  49,  50).  In  those  works,  however,  the  explicit  energy-gap 
functional  form  has  not  been  specified. 

Note  that  given  our  conclusion  about  the  n  energy-gap  symmetry,  we  have  in  MgB2  an 
interesting  configuration  where  the  symmetry  of  the  a  energy  gap  (s-wave)  and  that  of  the  n 
energy  gap  (i-wave)  are  different.  This  configuration  may  not  be  unique  to  MgB2.  We  note  that 
the  superconductors  2H-NbSe2  and  NbSe3  share  the  same  crystal  structure,  both  have  two  energy 
gaps,  and  comprise  alternating  Nb  and  Se  layers  analogous  to  the  B  and  Mg  layers  in  MgB2  [51], 
This  suggests  that  low  temperature  measurements  of  the  penetration  depth,  surface  resistance  and 
intennodulation  distortion  2H-NbSe2  and  NbSe3  would  be  of  great  interest. 
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5  Summary 

We  have  presented  a  phenomenological  model  for  calculating  the  energy  gaps  in  MgB2  at 
I=0K,  The  model  is  an  extension  of  the  BCS  theory  to  account  for  both  the  coulomb  repulsion 
and  the  considerable  anisotropy  in  the  elastic  properties.  The  calculations  support  our  previous 
hypothesis  that  the  n  energy  gap  is  nodal,  with  six  nodal  lines  £  =  6  (I  wave  in  spectroscopic 
notation),  in  contrast  to  the  common  assumption  that  the  n  gap  is  conventional  s  wave.  Our 
calculations  are  also  consistent  with  the  common  viewpoint  that  the  a  energy  gap  has 
conventional  s-wave  symmetry.  The  calculations  assume  a  Fermi  surface  that  comprises  of  two 
sheets  each  underlying  a  separate  energy  gap  [10],  and  no  coupling  between  the  phonons  and 
electron  seas  that  pertain  to  the  two  gaps.  All  the  required  input  parameters  for  the  model  are 
available  from  independent  data  and  calculations  in  the  literature.  Given  that  the  gap  equation  is 
highly  nonlinear  (Eqs.(4)  and  (B.10))  this  outcome  substantiates  the  self-consistency  of  the  input 
parameters  (Table  2). 

We  have  also  presented  a  new  analysis  of  recently  published  measurements  of  the  surface 
resistance  in  MgB2  [6],  This  data  clearly  exhibit  non-exponential  decay  at  low-temperatures,  a 
feature  inconsistent  with  an  s-wave  n  gap.  Our  analysis  shows  that,  in  fact,  the  measured  low- 
temperature  dependence  of  surface  resistance  follows  a  power  law,  as  would  be  expected  for  a 
nodal  energy  gap.  This  result  is  further  substantiated  by  a  survey  of  the  literature  on 
measurements  of  the  low-temperature  surface-impedance  variation  that  leads  to  the  conclusion 
that  there  is  no  unambiguous  evidence  for  the  exponential  dependence  as  predicted  by  an  s-wave 
symmetry. 

All  of  the  findings  here  support  strongly  our  previous  hypothesis,  based  on  our 
measurements,  that  the  n  energy  gap  is  nodal  with  £  =  6,  six  nodal  lines. 
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Appendix  A 

The  Phonon-Induced  Attractive  Interaction  VPH  (q) ,  Eq.(8) 

Equation  Section  (Next)The  parametrization  of  the  V PH  (q)  is  motivated  by  the 

corresponding  expression  for  an  isotropic  superconductor,  Eq.(6).  Consider  two  limits  of  the 
momentum  exchange,  i.e.,  in  the  _L  and  □  directions,  and  assume  that  in  (6)  the  only  direction- 

dependent  factor  in  the  amplitude  is  vSOUND .  For  these  special  instances  we  would  write 


r(q±)  = 

Consequently,  a  convenient  expression  that  applies  to  all  q  directions  and  coincides  with  the  two 
limits  in  Eq.(A.l)  is  the  interpolation 
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Since  the  n  -  gap  involves  vibrations  of  the  B  and  Mg  planes  against  each  other,  i.e.,  in  the  _L 
direction,  and  employing  (A.3)  we  write 
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In  (A.4)  we  used  the  substitutions  mentioned  after  Eq.(6).  By  the  same  token,  since  the  cr  -gap  is 
associated  with  in-plane  B  vibrations  in  the  □  direction,  we  arrive  at 


V,„(9)  =  - 
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In  Eqs.(A.4),(A.5),  the  symbols  Z?(_l_), /?([])  denote  the  bulk  moduli  in  the  _L  direction  and  the  □ 

plane,  respectively  (dimensionality  [»]  =  ££  3).  For  the  hexagonal-crystal-structure  the  elastic 
constants  are  given  by  [52,53] 

s(i ;□) — ,  s(±)=— - — - 

2  +  ccR  ccR  ( 2  +  ocR  ) 

C  +C  —2C 

aR= — - — ,  A  =  2(c1j1  +  C1j2)  +  4C1j3  aR+C33aR,  (A.6) 

^3,3  ^1,3 


where  the  two  moduli  in  Eq.(A.6)  obey  the  relation 


B{0)  B{±)  C,,(C,,+Cu)-2Cl, 
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Appendix  B 

Multipole  Expansions  and  the  Gap  Equation,  Eq.(4). 

Equation  Section  (Next) 

To  solve  the  gap  equation,  Eq.(4),  it  is  necessary  to  specify  the  band  structure  and  deal 
with  the  required  three-dimensional  integration.  These  matters  are  detailed  in  this  Appendix  to 
facilitate  reproduction  of  our  calculations  if  needed. 

Consider  first  the  band  structure.  Ab-initio  band-structure  calculations  yield  two  K 
Fermi-surface  sheets  (an  electron-like  and  a  hole-like),  adjacent  to  the  Brillouin  zone  facets  and 
of  complicated  topology,  and  two  nesting  a  tube-like  Fermi  surface  sheets  (hole  like)  aligned 
with  the  r-A  line  [9,54,  55].  In  our  calculations  we  consider  only  one  Fermi  sheet  for  each 
energy-gap  .For  the  K  Fermi-surface  sheet  we  approximate  the  hexagonal  Brillouin  zone  by  a 
cylinder  of  radius  and  half-height k/;z  =  nja  c  .  The  maximum  and  minimum  lateral 

2tt 

extensions  of  the  hexagonal  Brillouin-zone  geometry,  suggest  that  —j= - <  K 

where  a  a  is  the  lattice  constant  in  the  x  direction.  In  that  Brillouin  zone  cylinder,  the  n  Fermi- 

surface  sheet  extends  over  about  half  the  Brillouin-zone  height.  Since  we  average  the  gap 
equation  over  the  kz  variable,  see  below,  we  use  for  that  sheet  the  entire  Brillouin-zone  height. 
For  the  a  energy  gap,  the  chosen  Fermi  surface  is  a  single  tube  with  an  averaged  radius  K F  (a) 
[56].  These  considerations  suggest  the  following  approximate  band  energies 


■W< 


An 
3  a  n 


£x(K’kz)~Y^’  £a(K’kz)*Y^’  and  K<KF(n,a )  ~  —  <kz< 

TYl  ^  TYl  ^  ac 

and  ml,m*a,KF(n,  a)  are  effective  masses  of  the  K  and  cr  Fermi  surfaces, 
available  from  data  and  band-structure  calculations  [56,57]. 


— ,  (B.l) 

ar 


These  are 


Eq.(B.l)  implies  that  %(k)~  £(K)  and  that 
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K-K'\ 


=  KFy]  2(l-cos(p-p')) 


(B-2) 


where  (p,(p '  are  the  basal  azimuthal  angles  of  K,  K' ,  respectively,  and  KF  is  the  pertaining  basal 
Fenni  momentum  in  the  cylindrical  coordinate  system. 

To  carry  out  the  integrations  implied  by  (4)  the  multipole  expansions  of  the  interactions 
are  required.  These  are 


co 

V  PH  (  4  )  =  Z  V  f-H1’  (  ?  X  )  C0S  (  m  (<P-<P'))fm 

m  =  0 

co 

Vc  (?)  =  Z  Vcn>  (« ± )  cos (m (p - p ’))  fm 


m  =  0 


f  =  2-8 

8  m  n 


(B.3) 


where  m  denotes  the  multipole  order  and 
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(B.4) 


According  to  Eq.(l  1),  only  the  multipoles  m=0,3,6,  are  possible,. 


We  turn  now  to  the  RHS  of  the  gap  equation  (4),  which  we 


denote  here  by  F  {^k  j .  For  the 


approximate  band  energy  of  Eq.(B.  1)  we  have 
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\dk'F(${k  '),kz-kz>-p'): 


/  *  \  nla 

m  »  7 


J  dkz  d%'  J  d(p'  F (£;',kz-kz',(p - <p')  (B.5) 

-  V°c 


—  oo  0 


The  triple  integration  in  Eq.(B.5)  becomes  tractable  since  the  <p'  integrations  are  available 
analytically,  and  the  <p'  and  the  kz  '  variables  are  separable. 


Consider  first  the  angular  integration.  For  s-waves  ( m  =  0 )  the  (p  integration  in  is  trivial. 
For  m  >  0  the  angular  integration  gives  rise  to  the  integral 


2  K 

|  d(p 
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2  n 
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1  d 

1  ^  Jj 

(B-6) 


where  the  symbols  E,K  denote  the  complete  elliptic  integrals  of  the  first  and  second  kinds, 
respectively  [58,59]. 

The  integration  in  Eq.(4)  is  limited  by  the  step  functions  in  Eqs.(5),(l  1).  This 
integration  yields 


,(0) 


,  (m) 


U, 

f  ..  2 


U  +  yj  It  +  A 
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4f 


v  f  pp 

V  vf’y 


-K 


f  A2^ 


for  m  =  0 
for  m  =  3,  6.  (B.7) 


The  analytic  expression  for  G imi  involves  the  Meijer  G  functions  is  not  quoted  here.  Carrying 
the  f,  ’  and  (p  ’  integrations  results  in  that  the  <p,  f  dependencies  on  both  sides  of  the  gap 
equation,  Eq.(4),  cancel  out  exactly,  as  they  should. 

The  last  kz  integration,  Eq.  (1 1),  is  done  numerically.  For  the  two  symmetries  considered 
in  this  work  fr (kz)  —  \.  The  final  equations  for  the  two  unknowns  A  !'/:i  (T )  in  Eq.(ll)  are 
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obtained  by  averaging  both  sides  of  the  gap  equation  (4)  over  the  Brillouin  zone  kz  coordinate 
[36].  This  step  gives  rise  to  the  form  factors 


i  (m)  _ 

lPH,C 


7rjac  nfaf 


J  J  dkzdkz'V^c(kz-kz')  J  dkz. 


0  -nja 


(B-8) 


It  is  convenient  to  express  the  ensuing  equations  in  terms  of  dimensionless  unknowns. 
For  this  purpose  we  introduce  the  following  variables 


fid) D  p  = 


tlCQ, 


A(T  =  0/f)’ 
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(1.2) 
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(1,2) 


A(T  =  0K ) 


(B-9) 


where  A(r  =  0A")  is  the  pertaining  measured  energy-gap  {7T  or  cr).  Putting  all  these  pieces 
together,  the  coupled  nonlinear  equations  for  the  unknowns  A('n2)  are 
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(B.10) 

The  coupled  equations  (B.lO)are  solved  by  searching  for  a  solution  over  a  ‘reasonable’  range  of 
the  A(^  unknowns.  In  some  instances  there  are  no  solutions  except  for  the  trivial  one  A12;21  =  0  . 
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Table  1.  The  allowed  positive-parity  energy-gap  symmetries  for  hexagonal  crystal 
structure* 


Symmetry  notation 

r+ 

x  i 

r: 


Energy  gap 

A  (*) 

A„(r),  A„(r)|f, 

BZ 


BZ 


A(^;r)  =  A0(T)sin(6^) 

rsin(3^) 

[cos  (3^) 

*  The  symbols  KB z  and  kB/  denote  the  momentum  extensions  of  the  first  Brillouin  zone  in  the 


A(«,;7-)  =  A0(r)-^x. 

* BZ 


basal  plane  and  z  directions,  respectively.  These  arbitrary  nonnalizing  constants  are  introduced 
so  that  the  energy-gap  amplitude  A  0(T)  has  dimension  of  energy. 
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Table  2.  The  parameters  employed  in  Eq.(8) 


Parameter 

n  Energy  Gap 

<  Vf>/c 

0.00161 

z 

1/2 

n0  Vuc  (nm3) 

2 

B  (GPa) 

B(_L)=279.07 

m*/m0 

0.375 

A(T=0  K)(eV) 

0.002 

K  F(nm x) 

10.018 

a  Energy  Gap 

Reference 

0.00146 

60,54 

1 

See  text 

1 

See  text 

B(||)=587.78 

Eqs.  (12), (A. 6) 

0.9 

1.25.57 

0.0071 

1,2,37 

3.34 

1,2,37 

*  Vuc  =.0.25  a  a2a  CV3  is  the  primitive  unit-cell  volume,  where  we  used  a  a=0.3086  nm  ,  a  c  = 
0.3524  nm  [1],  Other  parameters  used  in  the  calculations  are:  ficoD  =750  kB  where  k  B  is  the 
Boltzmann  constant  [37,  55],  hojp  -7eV  ,  An=  0.002  eV,  A  a  =  0.007  eV  [1,9],  q  tf  =  0.7  (a  B)" 
'where  a  B  denotes  the  Bohr  radius[54],  ti c  =  197.3  eVxnm. 
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Figure  1  Normalized  IMD  vs  reduced  temperature  T/Tc  A  for  typical  YBCO[12],  and  •  for  MgE^ 
[3],  Solid  line  is  a  l/T2  dependence.  The  similarity  between  YBCO  and  MgB2  is  an  indication 
that  the  energy  gap  in  MgBi  contains  nodes. 
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Figure  2.  Main  plot:  A  A/A  vs.  temperature  for  MgB2  very  low  temperatures  [4],  The  inset  shows 
the  same  quantities  over  the  entire  temperature  range  <  Tc.  Several  samples  are  shown.  All 
demonstrate  the  increase  at  low  temperatures  attributed  to  Andreev  bound  states  [4],  This 
increase  is  only  visible  on  the  expanded  scale  of  the  main  plot. 
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1/T  (K'1) 

Figure  3.  Semilog  plot  of  Rs  =  R Measured  -  Rres  vs.  1  IT,  where  R Measured  denotes  the  measured 
effective  surface  resistance  of  the  350nm  film  by  the  JLab  group  [6]  and  Rres  is  defined  as  the 
temperature  independent  value  of  surface  resistance  at  low  temperatures.  The  solid  (red,  in  the 
online  version)  and  dashed  (black,  in  the  online  version)  are  plots  of  e  ' 1:1  where  A  are  accepted 
values  for  MgB2,  =  2.3  meV  and  Aa  =  7.1  meV,  respectively.  Neither  value  fits  the  data. 
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Figure  4.  Semilog  plot  of  measured  niobium  surface-resistance  vs  Tc  /T  from  JLab  [6].  •  is 

R Measured  values,  ■  R Measured  -Rres •  Solid  line  (red  in  online  version)  is  Eq.  (2)  with  the  accepted 
value  of  for  the  niobium  energy  gap,  A=  1.5  meV,  which  fits  the  data  over  two  orders  of 
magnitude. 
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Figure  5.  Rs  vs  T  of  the  JLAB  data  [6]  and  stripline  resonator  measurements  at  Lincoln 

Laboratory  of  MgEL  deposited  by  reactive  evaporation  [22].  All  data  are  scaled  by/  to  7.4  GFIz. 
▲  measurements  of  the  350-nm  thick  film  from  JLab,  ■  stripline-resonator  measurement  at 
Lincoln  Laboratory  of  film  deposited  by  reactive  evaporation  and  •  stripline  resonator 
measurements  at  Lincoln  Laboratory  of  MgEL  deposited  by  HPCVD.  The  lines  are  best  fits  to  the 
various  data  sets  with  the  expression  in  Eq.  (3)  with  An  =  2.3  meV,  or  explicitly 

Rs  =A(T /  Tc  +  .0026(7" ITcf)+B ,  where  A  and  B  are  taken  as  free  parameters.  The  values  of  A 

o 

and  B  are  as  follows:  for  the  long-dashed  line  A  =  4.04  and  B  =  3.87xl0'1 ;  for  the  solid  line  A  = 
1.30  and  B  =  27.3;  and  for  the  short-dashed  line  A  =  0.607  and  B  =  7.13  in  the  proper  units. 
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Figure  6.  Ratio  of  the  calculated  n  energy  gap  for  the  t  =  6  symmetry  A6  and  the  representative 
experimental  value  AK  =  2  meV  vs.  the  anisotropy  parameter  a  Eq.(9).  The  shaded  region  depicts 
the  range  of  experimentally  determined  anisotropy  parameters.  The  points  are  the  numerically 
calculated  values  and  the  line  is  a  spline  fit  to  the  points.  The  calculated  A6  is  within  an  order  of 
magnitude  of  the  experimental  value. 
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Figure  7  Ratio  of  the  calculated  n  energy  gap  for  the  t  =  0  symmetry,  Ao  and  the  experimental 
value  An  =  2  meV  vs.  the  anisotropy  parameter  a  Eq.  (9).  The  shaded  region  depicts  the  range  of 
experimentally  determined  anisotropy  parameters.  This  calculated  result  disagrees  with  the 
experimental  result  y  two  oders  of  magnitude.  Together  with  Figure  6  this  leads  to  the  conclusion 
that  the  i  =  6  symmetry  is  much  more  likely  for  the  n  gap.  The  points  are  the  numerically 
calculated  values  and  the  line  is  a  spline  fit  to  the  points. 
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Figure  8.  Ratio  of  the  calculated  a  energy  gap  for  the  t  =  0  symmetry,  A0  and  the  experimental 
value  Aa  =  7  meV  vs.  the  anisotropy  parameter  a  Eq.  (9).  The  shaded  region  depicts  the  range  of 
experimentally  determined  anisotropy  parameters.  There  is  no  nontrivial  solution  for  the  case  of 
l—  6  symmetry.  The  calculated  result  is  within  an  order  of  magnitude  of  the  experimental  value. 
This  indicates  that  the  a  gap  is  likely  a  conventional  s-wave  symmetry.  The  points  are  the 
numerically  calculated  values  and  the  line  is  a  spline  fit  to  the  points. 
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